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QCD with 2 flavours of massless colour-sextet quarks is studied as a possible 
walking-Technicolor candidate. We simulate the lattice version of this model at 
finite temperatures near to the chiral-symmetry restoration transition, to determine 
whether it is indeed a walking theory (QCD-like with a running coupling which 
evolves slowly over an appreciable range of length scales) or if it has an infrared 
fixed point, making it a conformal field theory. The lattice spacing at this transition 
is decreased towards zero by increasing the number iV) of lattice sites in the temporal 
direction. Our simulations are performed at Nf = 4,6, 8,12, on lattices with spatial 
extent much larger than the temporal extent. A range of small fermion masses 
is chosen to make predictions for the chiral (zero mass) limit. We find that the 
bare lattice coupling does decrease as the lattice spacing is decreased. However, it 
decreases more slowly than would be predicted by asymptotic freedom. We discuss 
whether this means that the coupling is approaching a finite value as lattice Nt is 
increased - the conformal option, or if the apparent disagreement with the scaling 
predicted by asymptotic freedom is because the lattice coupling is a poor expansion 
parameter, and the theory walks. Currently, evidence favours QCD with 2 colour- 
sextet quarks being a conformal field theory. Other potential sources of disagreement 
with the walking hypothesis are also discussed. 

We also report an estimate of the position of the deconfinement transition for 
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Nt = 12, needed for choosing parameters for zero-temperature simulations. 


INTRODUCTION 


The LHC at CERN is currently probing the Higgs sector of the Standard Model of high- 
energy physics. This sector is the least well understood part of the standard model, and the 
least satisfactory from a theoretical standpoint. Thus the study of extensions of the standard 
model with a more aesthetically compelling Higgs sector is timely. The observation of a light 
125 GeV) Higgs-like excitation at ATLAS and CMS, with properties consistent with the 
standard-mo del Higgs, puts constraints on any such model. We are especially interested in 
those models where the Higgs sector is strongly-coupled and the Higgs boson is composite. 

We are interested in QCD-like models - non-Abclian gauge theories with massless 
fermions and spontaneously broken chiral symmetry - where the pion-like Goldstone bosons 
play the role of the Higgs field, giving mass to the W ± and Z weak vector bosons through 
the Higgs mechanism. Here the Higgs boson is the remnant radial excitation. Such theories 
are called Technicolor models [1, 2 J]. Technicolor models, which are simply QCD scaled up so 
that f n ~ 246 GeV rather than f n & 93 MeV of regular QCD, are not phenomenologically 
viable. It has been suggested that Technicolor theories where the fermion content is such 
that the running gauge-coupling evolves very slowly over an appreciable range of mass scales, 
described as ‘walking’ rather than running, might be capable of overcoming such difficulties. 
Such theories are referred to as Walking Technicolor models . Because of their nature, 
the non-perturbative properties of such models are amenable to study using the simulation 
methods developed for Lattice QCD. It is such theories that we are interested in simulating. 

Candidate walking gauge theories typically have 2-loop /3-functions with a second, non¬ 
trivial, zero. If this behaviour remains true to all orders, this non-trivial fixed point controls 
the infrared properties of the theory, which is therefore a conformal field theory with a 
continuous spectrum. On the other hand, if the running coupling becomes so large that a 
chiral condensate forms before this would-be IR fixed point is reached, the theory is QCD- 
like with spontaneously broken chiral symmetry and Goldstone bosons separated by a mass 
gap from the rest of the spectrum. Because of the proximity of the would-be fixed point, a 
region where the coupling walks would be expected. 

We have concentrated our efforts on techni-QCD with 2 ffavours of massless Technicolor- 
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sextet techni-quarks. Since this is identical to QCD with 2 colour-sextet quarks scaled so 
that fa ~ 246 GeV, we will omit the prefix ‘techni’ from here on. This theory has been 
identified as a potential walking-Technicolor candidate (see for example jsl and references 
therein). It is asymptotically free and its 2-loop /3-function does have a non-trivial zero far 
enough from g 2 = 0 for the coupling to become sufficiently large for there to be a chance 
that chiral symmetry breaks before it is reached. If it is indeed QCD-like, it has 3 Goldstone 
bosons, the correct number to give masses to the Ws and Z, with none left over. In this 
sense, it is minimal. 

Other groups have studied/are studying this model using lattice techniques. The main 
contributors are Degrand, Shamir and Svetitsky 8-tl3] and the Lattice Higgs Collaboration 
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21], In addition we should mention some recent work by A. Hasenfratz and her collabora- 
22]. With the exception of one early paper by Degrand, Shamir and Svetitsky, 9] these 


tors 

have concentrated on the zero-temperature properties of this model. We study lattice QCD 
with 2 colour-sextet quarks at finite temperature. Our goal is to determine if the evolution 
of the coupling as lattice spacing a —> 0 is described by asymptotic freedom, and that chi¬ 
ral symmetry remains broken in this limit. Assuming that the chiral-symmetry-restoration 
transition is indeed a finite-temperature transition, increasing the temporal extent of the 
lattice, N t in lattice units, with the spatial extent N s » N t , and with temperature T fixed 
at the chiral phase transition temperature T x , takes a = l/N t T x towards zero. Thus the 
running coupling at this temperature, g x (a), should approach zero as Nt —> oo in a manner 
determined by the perturbative /3-function. If, on the otherhand, the theory is conformal, 
g x will approach a finite limit as Nt —> oo, characterising a bulk transition. Similar argu¬ 
ments should apply to the deconfinement transition at gd- However, as we have determined, 
deconfinement occurs at a much stronger coupling than chiral-symmetry restoration. For 
the N t values we have considered (N t < 12), g c j is too large for its evolution to be controlled 
by the perturbative /3-function. 

Recent use of this method, to search for the lower bound (in Nf) of the conformal window 
for QCD with many fundamental quarks 23 26] have shown it to complement step-scaling 
methods. The reader should consult the references in these papers for the history of such 
studies. Although these methods are, in principle, straight forward, these recent papers 
indicate that they are not so easy to implement in practice. 

Our earlier studies at AA = 4, 6, 8 { 27 , 28] were consistent with the evolution of g x between 
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N t — 6 and N t — 8 being described by the 2-loop /3-function. We have extended our 
simulations to N t = 12. In addition, we have covered the neighbourhood of the chiral 
transition for N t — 6 and N t — 8 with more closely spaced values of /3 = Q/g 2 to determine g x 
more precisely. In addition, we have determined the position of the deconfinement transition 
for N t = 12 for one mass value. Preliminary versions of the results presented in this paper 
have been presented at lattice conferences 29j, [30]. 

While the observed change in f3 x = 6 jg 2 x between N t — 6 and N t — 8 is consistent 
with that predicted using the 2-loop /3-function, that between N t — 8 and N t = 12 is only 
about half the predicted value. At face value, this suggests that j3 x could be approaching 
a finite limit, which would mean that this theory is conformal. However, we need to be 



this and other potential sources of systematic errors in our approach, later in this paper. 
Our results make it important to perform further studies of lattice QCD with 3 massless, 
colour-sextet quarks, to determine if it approaches its expected asymptotic behaviour at 
N t = 12, which would be qualitatively different from that observed for 2 flavours. 

In section 2 we discuss our methods of simulation and analysis. Section 3 gives the 
results of our simulations at N t — 6, 8 and 12 near the chiral transition. We compare these 
results with perturbative predictions in section 4, and discuss improvements. In section 5 
we analyse simulations at N t = 12 in the neighbourhood of the deconfinement transition. In 
section 6 we discuss our results, try to draw conclusions, and indicate directions for future 
studies. 


II. METHODS 


This section largely repeats discussions given in earlier publications [27|, [28], and is in¬ 
cluded here for completeness. We use the simple Wilson plaquette action for the gauge 
fields: 


s 9 = PE 


r Re(Tr UUUU) 


( 1 ) 
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Here the gauge fields U on the links are in the fundamental representation of SU (3) co i our . 
We use the unimproved staggered-fermion action for the quarks: 

'Nf/4 


S, = £ 


sites 


Y Y f [ft+rn\ft f 
f =i 


( 2 ) 


where ft = E,, with 


D^ix) = + ft) - U^ftx - jl)ft(x - £)], 


( 3 ) 


where [A 6 ) is the sextet representation of U. i.e. the symmetric part of the tensor product 
U ®U. Reasons for this choice have been discussed in our earlier publications. When Nf is 
not a multiple of 4 we use the fermion action: 


Sf=Y xHlP+ m }[-ft+m]} Nf/8 X- 


( 4 ) 


sites 


The operator which is raised to a fractional power is positive definite and we choose its 
positive-definite root. This yields a well-defined operator. We use the RHMC method 


for our simulations 


341 . where the required powers of the quadratic Dirac operator are 


replaced by diagonal rational approximations, to the desired precision. By applying a global 
Metropolis accept/reject step at the end of each trajectory, errors due to the discretization 
of molecular-dynamics ‘time’ are removed. 

The canonical partition function for a field theory at finite temperature T is realized 
by evaluating the functional integral for Euclidean time where the time is restricted to 
an interval 1/T with periodic boundary conditions on the boson fields and antiperiodic 
boundary conditions on the fermion fields. Space is kept infinite. On a lattice of lattice 
spacing a, this means using a lattice of temporal extent N t in lattice units where N t a = 
1/T. The spatial extent of the lattice N s » N t . For lattice QCD with sextet quarks, 
if deconfinement and chiral-symmetry restoration are finite-temperature transitions, then 
the associated temperatures at which they occur, T d and T x respectively, should be fixed, 
independent of a, for a small enough. Thus measuring the couplings at either of these 
transitions as N t is varied, gives g(a) for a sequence of as which approaches zero as Nt —» oo. 
As it turns out /3d and hence g d lies in the strong-coupling domain, outside the regime where 
perturbation theory is likely to be valid, for N t = 4, 6, 8,12 and any other N t which we 
are likely to consider in the near future. We therefore concentrate our efforts on the chiral 
transition, which occurs at much weaker couplings. If, on the other hand, QCD with 2 
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colour-sextet quarks is conformal, the chiral transition would be a bulk transition. In this 
case g x would approach a non-zero limit for large N t , and the whole region of broken chiral 
symmetry would be a lattice artifact, disconnected from the conformal held theory at weaker 
coupling. 

If QCD with 2 colour-sextet quarks is QCD-like, the approach of g x to zero is described 
by asymptotic freedom expressed in terms of the /3-function. Through 2 loops this is given 
by: 

13(g) = -big 3 - b 2 g 5 . (5) 

Expressed in terms of (3 = 6/g 2 , the evolution of the coupling when the lattice spacing is 
scaled by A is given by 

A 13(13) = /3(a) - (3(\a) = (12^ + 72 b 2 //3) ln(A) + 0( 1//3 2 ), (6) 


where for Nf flavours of colour-sextet quarks: 

h = (H - y N >) ^ 16,r2 

b 2 = (l02 - ™N,) /(16 7T 2 ) 2 . (7) 


The chiral transition occurs at that value of (3 (f3 x ) at which the chiral symmetry is 
restored and beyond which the chiral condensate ('tpij)) vanishes, for massless quarks. Of 
course, in lattice simulations, we need to run at (small but) finite mass, and extrapolate 
to zero quark mass. It is not, however, practical to run at masses small enough for the 
condensate to accurately determine (3 X directly. We therefore estimate (3 X from the peaks 
in the chiral susceptibilities, or rather in the disconnected part of the chiral susceptibilities. 
This is given by: 


Xyy = V (WO 2 ) - WO 2 


( 8 ) 


where the () indicates an average over the ensemble of gauge configurations and V is the 
space-time volume of the lattice. Since we use stochastic estimators for 'W, we need at 
least 2 estimators per configuration. The Erst term must include only contributions which 
are off-diagonal in the noise, to obtain an unbiased estimator. We, in fact, use 5 stochastic 
estimators at the end of each trajectory giving 10 estimates for x^ P er configuration. 



III. SIMULATIONS OF QCD WITH 2 FLAVOURS OF COLOUR-SEXTET 

QUARKS AT N f = 6,8,12 
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Here we describe our simulations with 2 colour-sextet quarks on lattices with Nf = 6,8 
and 12. For Nf = 6 and 8 we have extended the simulations of our earlier papers, where 
we simulated at f3 spacings of 0.1, through the chiral transition region. For the lowest mass 
(m = 0.005) at Nf = 6, we have covered the vicinity of the chiral transition at f3 spacings of 
0.02, and with increased statistics. At N t — 8 we have also covered the vicinity of the chiral 
transition with (3 spacings of 0.02 for all masses, including a new smaller mass (m = 0.0025). 
We have performed new high-statistics simulations at N t = 12, covering the region of the 
chiral transition with (3 spacings of 0.02, for all masses. Preliminary results of these new 
simulations have been presented at Lattice 2011, Lattice 2012, Lattice 2013 and Lattice 
2014. In each case, we simulate using the RHMC algorithm with trajectory length 1. Most 
of our simulations have been performed on lattices with N s = 2 N t . 

A. N t = 6 

Our simulations at N t = 6 are performed on 12 3 x 6 lattices with quark masses m = 0.005, 
m = 0.01 and m = 0.02. To more accurately pinpoint the peak of the chiral susceptibility 
at the lowest mass (m = 0.005) we have covered the region in the neighbourhood of the 
chiral transition 6.5 < /3 < 6.7 at [3 spacings of 0.02. At each of these /3s we have performed 
runs of 100,000 trajectories. Outside of this interval, and that near the deconfinement 
regime, we employ [3 spacings of 0.1 and 10,000 trajectories per f3. The chiral susceptibilities 
(equation [8]) from these new runs and those for m = 0.01 from our earlier work are plotted 
in figure [D 

Here we attempted to determine the position of the peak of the m = 0.005 susceptibility 
using Ferrenberg-Swendsen interpolation, but were unable to obtain consistent results. We 
therefore chose to fit the ‘data’ with a smooth curve: 

Xw = a-b(/3- /3 X ) 2 -c(/3- f3 x ) 3 . (9) 

The rational for this simple form is so that we can use the same form for each N t . The second 
term is to give a simple parabolic fit to the peak. The third term is necessary because, as is 
obvious for the larger N t s, the susceptibility is not symmetric around the peak. The value 


12 3 x 6 lattice 
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FIG. 1: Chiral susceptibilities on a 12 3 x 6 lattice, with Nf = 2. 

obtained for /3 X from a fit using all /3s in the range 6.4 < (3 < 6.7 is j3 x = 6.611(3) for a fit 
with \ 2 /d.o.f = 1.55, which is acceptable. This fit is shown, superimposed on the data in 
figure O 


B. N t = 8 


We have extended our simulations at N t — 8 on a 16 3 x 8 lattice. For the 3 masses 
considered in our earlier work, m = 0.005, m = 0.01 and m = 0.02, we have increased 
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12 3 x6 lattice N f =2 m=0.005 



P 


FIG. 2: Chiral susceptibilities on a 12 3 x 6 lattice, with Nf = 2 m = 0.005. The curve is the fit 
described in the text, with a = 26.0489, b = 407.578, c = 978.166 and (i x = 6.61143. 

the number of values in the neighbourhood of the chiral transition, 6.6 < (3 < 6.8, by 
simulating at /3s separated by 0.02 compared with our previous 0.1. We have increased our 
statistics to 50,000 trajectories at each (/3, m) in this range. In addition, we have simulated 
at a new, lower mass, m = 0.0025. Again we have covered the range 6.6 < (3 < 6.8 with 
/3s spaced by 0.02, performing runs of 100,000 trajectories at each f3. Outside this range we 
performed a run of 20,000 trajectories at f3 = 6.5, and runs of 10,000 trajectories at (3 = 6.9 
and (3 = 7.0, for this smallest mass. Figure [3] shows the chiral susceptibilities for these runs. 
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16 3 x8 lattice, N f = 2 
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FIG. 3: Chiral susceptibilities on a 16 3 x 8 lattice with Nf = 2. 

To estimate the position of the peak of the chiral susceptibility for m = 0.0025, we 
first consider using Ferrenberg-Swendsen interpolation of the chiral susceptibilities. This is 
possible, since the distributions of plaquette values for adjacent /3s in the range 6.6 < f3 < 6.8 
show significant overlap. Here we performed extrapolations from the susceptibilities for 
(3 = 6.66, f3 = 6.68, (3 = 6.70, (3 = 6.72 and f3 = 6.74, and looked for consistency in our 
predictions. The best consistency we found was between extrapolations from (3 = 6.68, which 
predicted a peak at (3 = 6.691(24), and f3 = 6.70, which predicted a peak at (3 = 6.689(5). 
Combining these we obtain a prediction (3 X = 6.69(1). 




11 


16 3 x 8 lattice N f =2 m = 0.0025 
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FIG. 4: Chiral susceptibilities on a 16 3 x 8 lattice with Nf = 2 and m = 0.0025. The curve is the 
fit described in the text, with a = 34.6359, b = 940.687, c = 2716.49 and /3 X = 6.70613. 

A second estimate comes from fitting our susceptibilities to the form we used for N t = 
6 (equation [9]). For m = 0.0025, fitting to this cubic polynomial over all points in the 
range 6.5 < /3 < 6.8 yields /3 X = 6.706(1) for the value of /3 at the peak. This fit has 
X 2 /d.o.f = 0.55, which is excellent. Figure @] shows this fit superimposed on the ‘data’. 
Performing a similar fit to the susceptibilities at m = 0.005 over the range 6.5 < j3 < 6.9 gives 
f3 x = 6.701(4) with y 2 fd.o.f = 0.85. Fitting the susceptibilities for m = 0.01 over the range 
6.5 < (3 < 6.9 yields (3 X = 6.693(8) with y 2 fd.o.f = 1.36 while fits to the susceptibilities 
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for m = 0.02 over the same range predicts [3 X = 6.71(1) with x 2 jd.o.f = 0.29. A word of 
caution is due concerning the fits for the 2 largest masses. In both these cases, the measured 
susceptibility is statistically flat over an appreciable neighbourhood of the transition. The 
positions of these 2 peaks is thus determined largely by the outlying points on the fits, and 
should therefore not be taken too seriously. The main reason for performing these high-mass 
fits is to get an estimate of the height of these peaks from the parameter a in the fits. 


X 


6 


~A- 

X 


16 3 x 8 lattice 



FIG. 5: Peak of chiral susceptibility as a function of mass, with fit to critical scaling form Xmax = 
on a 16 3 x 8 lattice. 


For a second-order phase transition, the value Xm at the peak, Xmax is expected to scale 
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with mass as: 

Xmax = Am 115 - 1 . (10) 

If the chiral transition is a finite-temperature transition, it is expected to lie in the 0(2) or 
0(4) universality class where the critical exponent 6 ps 4.8. If it is a bulk transition, which 
is expected to be first order, 5 = oo. The best fit to equation flOl gives 5 = 4.1(1) and has 
X 2 /d.o.f. = 9. Figure [5] shows this fit superimposed on the values of Xmax taken from the 
values of a in the susceptibility fits. Clearly the reason that the estimated quality of the 
fit (reduced x 2 ) is poor is because the systematic errors associated with choosing a as an 
estimate for Xmax have been ignored, whereas, especially for the larger masses, these clearly 
dominate. 

At N t = 8, we have also performed simulations with m = 0.0025 on a 24 3 x 8 lattice at 
3 (3 values, near and above the chiral transition, to check for finite volume effects. Values 
of the chiral susceptibilities, which should be most susceptible to finite volume effects, are 
given here, along with their values on a 16 3 x 8 lattice in square brackets. For [3 = 6.7, 
XXi> = 34.4(4) [35.0(5)], for (3 = 6.76, Xyy = 30.9(4) [31.5(4)], while for /3 = 6.9, xh> = 
7.0(2) [7.0(2)]. These results are in good-enough agreement for us to conclude that finite 
volume effects are small at the masses we use. 

C. N t = 12 

We perform simulations on a 24 3 x 12 lattice at masses m = 0.0025, m = 0.005 and 
m = 0.01, in the neighbourhood of the chiral transition. At the smallest mass, m = 0.0025, 
we perform simulations at f3 values spaced by 0.02 over the range 6.6 < (3 < 6.9. For 
6.6 < (3 < 6.66 we perform runs of 50,000 trajectories at each (3. For the range 6.68 < (3 < 6.9 
we perform runs of 100,000 trajectories per (3. At (3 = 6.5, we run for 25,000 trajectories, 
while for f3 = 7.0, (3 = 7.1 and f3 = 7.2, we perform runs of 10,000 trajectories per (3. At 
mass m = 0.005 we again cover the interval 6.6 < (3 < 6.9 at increments of 0.02. For f3 = 6.6, 
we run for 100,000 trajectories. For the rest of the interval i.e. 6.62 < (3 < 6.9 we run for 
50,000 trajectories per [3. At (3 — 6.5 we run for 25,000 trajectories, while for [3 = 6.4, 
/3 = 7.0, (3 = 7.1 and f3 = 7.2, we run for 10,000 trajectories per f3. At m = 0.01 we perform 
simulations of 25,000 trajectories per f3 at /3s spaced by 0.02 over the range 6.6 < (3 < 6.9. 
At (3 — 6.3, 6.4, 6.5 we run for 25,000 trajectories per f3, while for [3 = 6.2 we run for 12,500 
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trajectories. For (3 = 6.0, 6.1 and for (3 = 7.0, 7.1, 7.2 we perform runs of 10,000 trajectories. 
We also run for 50,000 trajectories per (3 for (3s spaced by 0.02 over the range 5.7 < (3 < 5.9, 
which is in the neighbourhood of the deconfinement transition. These runs at (3 < 6 will be 
discussed further in section 5. 

Figure [6^ shows the chiral condensates (( , 0'0)), as functions of f3 for all 3 masses m = 
0.0025, 77i = 0.005, m = 0.01. These are bare (lattice) quantities. However, at non-zero 
mass, if we expand in powers of the quark mass m, the coefficient of m in physical units 
diverges as 1/a 2 as a —* 0, and should be regularized. We therefore subtract part of this 
divergence using the prescription adopted by the Lattice Higgs Collaboration, where the 
subtracted chiral condensate is defined by: 


{^)sub = (V’V’) 


d 


m v 


dmy 


-(ipip) 


1 


my=m 


(ii) 


where my is the valence-quark mass. What we observe is, that while the unsubtracted chiral 
condensate shows indications that it will vanish in the chiral (m —> 0) limit for /3 sufficiently 
large, the subtracted chiral condensate shows this vanishing more clearly. However, even 
the subtracted chiral condensate does not yield an estimate for (3 X which is accurate enough 
for our purposes. We thus turn to using the peaks of the chiral susceptibilities, extrapolated 
to zero mass as our estimates for f3 x . 

Figure [7] shows the chiral susceptibilities, defined in equation |8] extracted from our mea¬ 
surements of (t/jt/j) (5 per trajectory) in our simulations on 24 3 x 12 lattices for masses 
m = 0.0025, m = 0.005 and m = 0.01. Here, the distributions of plaquette values for adja¬ 
cent (3s have insufficient overlap to even attempt using Ferrenberg-Swendsen interpolation 
to estimate the positions of the peaks in the susceptibilities for the 3 masses. The suscep¬ 
tibilities for Til = 0.0025 show a clear peak, those for m = 0.005 show some indication of 
a rather flat peak, while those for m = 0.01 show little evidence for any peak. In order to 
extract an estimate of the positions of these peaks, we use the fitting form used for N t — 6 
and 8 (equation |9|), which makes maximal use of the ‘data’. The fit to the m = 0.0025 
susceptibilities for all points in the range 6.5 < (3 < 6.9 yields (3 X = 6.768(2) in an accept¬ 
able fit with x 2 /d-o.f. = 1.68. This fit is superimposed over the measured susceptibilities 
in figure El A fit to the m = 0.005 ‘data’ over the same interval predicts [3 X = 6.745(6) 
with x 2 /d-o.f. = 0.91, while a fit to the tti = 0.01 data also over the same range gives 
(3 X = 6.70(3) with y 2 /d.o.f. = 1.52. These last 2 fits should not be considered too seriously, 
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FIG. 6: a) Unsubtracted chiral condensates (tpip) as functions of (3 for masses m = 
0.0025, 0.005,0.01. b) Chiral condensates subtracted using the Lattice Higgs Collaboration’s 

prescription, as functions of /3 for masses m = 0.0025,0.005,0.01. 
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24 3 x 12 lattice, N f = 2 
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FIG. 7: Chiral susceptibilities on a 24 3 x 12 lattice with Nf = 2. 

because these peaks are defined by outlying rather than central points, owing to the flat¬ 
ness of the distributions. Their main purpose is to yield an estimate for the values of the 
parameter a at their peaks. 

We also consider the scaling properties of the susceptibility peak with mass using the 
scaling form of equation [10] With only 3 points, there is only one degree of freedom. Our 
fit gives the critical exponent 5 = 2.98(4) with \ 2 jd.o.f. = 0.29. This is in agreement with 
the mean-field (free-field) critical exponent 5 = 3 , rather than the expected 0(2) or 0(4) 
critical exponent 5 ~ 4.8 or first-order scaling 5 = oo. The graph of our data - the values 
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24 3 x 12 lattice N f = 2 m = 0.0025 





FIG. 8: Chiral susceptibilities on a 24 3 x 12 lattice with Nf = 2 and m = 0.0025. The curve is the 
fit described in the text, with a = 25.691, b = 440.589, c = 1082.3 and f3 x = 6.76801. 

of a from our fits - with this fit superimposed is given in figure [UJ 

IV. COMPARISON WITH PERTURBATIVE PREDICTIONS 

Here for consistency we only consider the values of (3 X obtained from fitting our chiral 
susceptibilities to the form given in equation [9] The errors given in the previous section do 
not include any systematic errors due to this rather arbitrary choice of a fitting function or 
to the selection of the range of (3 values over which the fits are performed. We will assume 
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24 3 xl2 lattice 



FIG. 9: Peak of chiral susceptibility as a function of mass, with fit to critical scaling form Xmax = 
Am x ^~ x on a 24 3 x 12 lattice. 

0.01 as a conservative estimate of these systematic errors, at least for the lightest masses. 
Because our measurements are inadequate to reliably determine if there is any significant 
mass dependence on the position of the peaks, we take our measurement of the position of 
the peak for the lightest mass as our estimate for the position of the transition in the chiral 
limit. We assume that our estimate of the systematic errors is large enough to encompass 
the shift in the positions of the peaks as we approach the chiral limit. Table Hj summarises 
the results of our determinations of the positions of the chiral transitions (/ 3 X ) from this and 
previous publications, as well as those of the positions of the deconfinement transitions ((3d) 
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from previous publications and the next section. 


N t 

Pd 

Px 

4 

5.40(1) 

6.3(1) 

6 

5.54(1) 

6.61(1) 

8 

5.65(1) 

6.71(1) 

12 

5.81(1) 

6.77(1) 


TABLE I: Nf = 2 deconfinement and chiral transitions for Nt = 4, 6,8,12. 


Equation [6] gives the perturbative asymptotic freedom predictions for the changes in /3 X 
between different N t s. In particular: 

P x (N t = 8)~P x (N t = 6) » 0.087 

P x {N t = 12) - p x {N t = 8) » 0.122 (12) 

P x {N t = 12) - (3 x (N t = 6) » 0.209 . 

The reason for including the third equation is because, if this is a finite temperature tran¬ 
sition, the fermion mass in physical units for m = 0.005 at N t — 6 and for m = 0.0025 at 
N t = 12 are identical, so that one might hope that any error due to not extrapolating to the 
chiral limit might be minimized. Our measurements (from table H]) give 


P x (N t = 8) - P x (N t = 6) = 0.10(1) 

P x (N t = 12) - p x (N t = 8) = 0.06(1) (13) 

Px(Nt = 12) — P x (Nt = 6) = 0.16(1). 


If taken at face value, these favour the conformal option, where the fact that P x {N t = 
12) — ft x (N t = 8) is roughly half the value predicted by asymptotic freedom could indicate 
that P x is approaching a non-zero limit. 

We need to be cautious, since the lattice coupling and hence 1/P is known to be a poor 
expansion parameter. That is, higher order terms for the expansion of any quantity in 
powers of g 2 tend to be large. The simplest improved choice of /3s is the tadpole-improved 
P of Lepage and Mackenzie 31]: 


P = -(Tr a UUUU)P . 

o 


( 14 ) 
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For connection with Lepage-Mackenzie, /3 = Q/g 2 = 6/47ro;. (Note that for staggered 
fermions, tadpole improvement of the fermion determinant is equivalent to rescaling the 
fermion mass and can thus be ignored, since we are interested in the limit m —> 0.) The 
plaquette in the above equation should be evaluated at /3, on a lattice which is at zero 
temperature. Since T = 0, in practice means on an Nf lattice for which /3 « f3 ( j(N t ), this 
would require simulating on lattices much larger than any we contemplate. For this reason, 
we use the finite temperature plaquettes from our simulations in this equation. This yields 
P x (N t = 6) = 4.48(1), P x (N t = 8) = 4.58(1) and P x (N t = 12) = 4.65(1). This gives 

P x (N t = 8)-P x (N t = 6) = 0 . 10 ( 1 ) 

Px(N t = 12) - p x (N t = 8) = 0.07(1) (15) 

P x {Nt = 12) — P x (N t = 6) = 0.17(1). 

compared with the perturbative prediction: 

P x (N t = 8) - p x (N t = 6) « 0.083 

P x (N t = 12) - 0 x (N t = 8) « 0.117 (16) 

P x (N t = 12) - p x (N t = 6 ) « 0.200 . 

While this is an overall improvement, it is insufficient. Choosing instead f3y, the (3 associated 
with the inter-quark potential, makes little difference. Here f3y = 6 /gy = Q/Anay. We use 
the relation between a and ay from Lepage-Mackenzie to obtain (3y. The reason that going 
from (3 to (3y does not make much difference is because, as noted by Lepage-Mackenzie, the 
perturbative relation: 

P = (3y + 2.245 + 0(1/p v ) (17) 

(where we have chosen the momentum scale at which we measure Py to be n/a), is a good 
approximation. With such a constant shift, differences in /3s are left unchanged. In the 
prediction, based on the perturbative /3-function, for the /3s (Pys) we consider, the 2-loop 
contribution is small, so the perturbative predictions for differences in /3s and /3ys are almost 
the same. Here we see that replacing the lattice coupling with an improved coupling such 
as gy does not significantly affect changes in /3, since g is small enough that the 2-loop 
contribution to the /3-function is significantly smaller then the 1-loop contribution, for both 
the original and improved schemes. Remember that the 1- and 2-loop coefficients in the 
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/3-function are scheme independent. Hence changing from lattice to improved couplings 
will not significantly improve agreement between measured and predicted running of the 
couplings. Of course, choosing a much smaller momentum scale for /3y, driving it towards 
the perturbative fixed point, could improve agreement for /3 x (N t = 12) — /3 x (N t = 8), but 
would be difficult to justify. 


However, it is well-known that even with tadpole improvemen 


bation theory for staggered fermions is still badly behaved 
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of the gauge links, pertur- 
331. The reason is another 


form of tadpole, the ‘doubler tadpole’ due to flavour(‘taste’)-mixing responsible for taste 
breaking. Unfortunately for us, before an improved perturbation theory could be developed 
for staggered fermions, interest shifted to improved staggered fermions designed to minimize 
taste breaking, making perturbation theory better behaved. 

It is interesting to note that the lack of significant improvement using /3 (or /3y or 
instead of /3 has also been noticed by 23] in their studies of QCD with 8 fundamental 
quarks using an improved staggered-quark action, so perhaps it is a property of theories 
with slowly varying running couplings and not an artifact of using unimproved actions. This 
contrasts with the one system where there are precise measurements of the finite temperature 
transition for a large range of N t values, namely pure SU( 3) Yang-Mills theory (quenched 
QCD) transcribed to the lattice using the Wilson (plaquette) action. For this system using 
an improved coupling greatly improves the agreement between the measured N t dependence 
of the critical coupling /3 C and the prediction from the 2-loop /3-function, as shown in table [Til 
For these calculations we used the values of / 3 c (N t ) for quenched lattice QCD given in the 


N t N' t 

\J3cW)-Pc{Nt)\lattice 

[hW)-Pc{N t )\lattice 

\MN;,)-MNt)\2-ioo P 

\/3 c (Ni)-/3 c (N t )]i!-i oor> 

6 8 

34% 

19% 

8 12 

24% 

13% 

12 18 

17% 

8% 


TABLE II: Difference between the changes in the lattice /3 C with change in Nt and the prediction 
from the 2-loop perturbative /3-function, compared with the same quantity using the improved /3c, 
for quenched lattice QCD. 


recent publication 35]. 
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V. THE N t = 12 DECONFINEMENT TRANSITION. 

In section 3 we mentioned that we have extended our simulations on a 24 3 x 12 lattice 
at m — 0.01 into the neighbourhood of the deconfinement transition. Although fid is too 
small for its evolution to be governed by perturbation theory, knowledge of its value as a 
function of N t is necessary when choosing fi values for zero temperature simulations. It 
is also useful to know the value of the deconfinement temperature as well as the chiral- 
symmetry restoration temperature in physical units i.e. in terms of f n ~ 246 GeV. We 
chose to simulate this regime at only one quark mass so that we could devote most of our 
resources to the proximity of the chiral transition. For this same reason we chose the highest 
of our 3 masses. Here we are relying on the observation based on our studies at smaller N t s 
that fid depends only weakly on the quark mass. 

The position of the deconfinement transition is determined by the point below which 
the Wilson Line (Polyakov Loop) becomes very small. Figure flOl is a plot of Wilson Lines 
against fi for our 24 3 x 12 simulations at all 3 masses. For the m = 0.01 plot, we notice 
that the Wilson Line is near zero for small fis and then jumps to a value appreciably greater 
than zero at fi ~ 5.8. 

To determine the position of this deconfinement transition more precisely, we ran for 
50,000 trajectories per fi for fi values spaced by 0.02 over the interval 5.7 < fi < 5.9. 
Figure |TT] shows histograms of the distribution of magnitudes of the Wilson Line, for fi 
values close to the deconfinement transition. Note the qualitative difference between these 
histograms for fi < 5.80 and those for fi > 5.82. Those histograms for the lower range of fis 
are peaked at ~ 0.01, while those for fis in the upper range are peaked at ~ 0.03. From this 
we estimate that the deconfinement transition occurs at fi = fid = 5.81(1). The transition is 
very abrupt, suggestive of a first-order phase transition. This value for fi d has been entered 
in tabled] in section 4. 


VI. DISCUSSION AND CONCLUSIONS 

We use studies of scaling of the assumed Unite-temperature chiral transition of lattice 
QCD with 2 colour-sextet quarks to attempt to determine whether this theory is conformal 
or walking. This provides an alternative to the use of step-scaling methods to achieve this 
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24 3 xl2 lattice 



P 


FIG. 10: Wilson Lines as functions of (3 on a 24 3 x 12 lattice. These are traces of products gauge 
links in the colour-triplet representation of SU(3) co i our . 


goal. Recent extensive step-scaling studies using improved staggered quarks [21] indicate 
that this theory walks, while similar studies using improved Wilson quarks [22] present 
evidence for a fixed point, which would mean that the theory is conformal. Hence use of a 
different method to try and determine which is the correct behaviour is warranted. 

We simulate lattice QCD with 2 flavours of colour-sextet quarks on lattices with N t — 6,8 
and 12 in the vicinity of the chiral-symmetry restoration transition, to accurately determine 
the value (3 X of f3 = Q/g 2 at that transition. Under the assumption that this is a finite- 
temperature phase transition, which thus occurs at a fixed temperature in physical units, 
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24 3 xl2 lattice m = 0.01 


m 

H 

H 

> 

H 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 

|WILSON LINE| 


FIG. 11: Histograms of magnitudes of the Wilson Line for /3 values close to the deconfinement 
transition on a 24 3 x 12 lattice with m = 0.01. 

this gives the running of the gauge coupling g as the lattice spacing a is decreased. We 
compare the evolution of the coupling g x at this transition with the prediction from the 2- 
loop perturbative /3-function. What we find is that the change in j3 x = Q/g x is in approximate 
agreement with the perturbative prediction between N t — 6 and 8, but is smaller by about 
a factor of 2 than the prediction between N t = 8 and 12. This suggests that this chiral 
transition is a bulk transition for which (3 X would approach a finite value in the large N t limit. 
In this case this theory would be a conformal field theory, and not the desired walking theory. 
On the other hand, the scaling of the susceptibility peaks with mass suggests that the chiral 
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transition is second order, whereas the simplest scenario for a bulk transition suggests that 
it should be first order. However, the fact that the critical exponent 5 ~ 3 suggests that the 
transition is mean-held, which would be more likely for a true 4-dimensional and hence bulk 
transition, than for a finite-temperature and hence quasi-3-dimensional transition, which 
should be in the universality class of the 3-dimensional 0(2) or 0(4) spin model with 6 & 4.8. 

In our initial comparison, we used the prediction in terms of the bare lattice coupling 
g(a), which is known to be a poor expansion parameter. We therefore looked at the tadpole- 
improved coupling g, which is supposed to be a better expansion parameter, as well as the 
related coupling gy which is related to the heavy-quark potential. Both these improvements 
make a slight improvement in the running of the coupling, but not nearly enough to produce 
agreement with the perturbative results. (Based on our experience with quenched QCD we 
would have expected an agreement to within 20% or probably better.) We do note, however, 
that even when such improvements are applied, perturbation theory with unimproved stag¬ 
gered quarks does not work well. In addition, the examples where use of such couplings ( g , 
gy or gjfs) has improved the behaviour of lattice perturbation theory have been for theories 
such as QCD with quarks in the fundamental representation of the colour group, where the 
massless theory has only one mass or length scale, the scale associated with confinement 
and chiral symmetry breaking. For QCD with colour-sextet quarks, we have shown that 
the scales of confinement and chiral-symmetry breaking are very different. It is unclear if 
methods which work for a single-scale theory will continue to work for a two-scale theory. 

We note that use of staggered quarks and, in particular, unimproved staggered quarks, has 
potential difficulties because flavour symmetry (sometimes referred to as taste symmetry) 
is explicitly broken, and is only restored in the continuum limit. This means that in the 
chiral limit at non-zero lattice spacing, there are less massless degrees of freedom than in the 
continuum limit. Because of this, it is possible that a theory with an infrared fixed point and 
hence conformal, could appear to be walking. Hence, care needs to be taken in taking the 
continuum and chiral limits. Using an improved action can reduce these problems. However, 
improved staggered fermions have their own difficulties. Improving the staggered-fermion 


action can introduce extra phases which are lattice artifacts. This was_o 
the apparent bulk transition in QCD with 12 fundamental quarks 


36 


37 


Dserved in studies of 
and provides extra 


complications for studying the finite temperature behaviour of QCD with many fundamental 
quarks 
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26)] . Improvement is designed to produce actions whose weak coupling physics 
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is closer to that of the continuum theory. However, it does not necessarily produce better 
behaviour in the intermediate or strong-coupling domain. 

Another word of caution is necessary concerning our results. We have only used one 
spatial size (24 3 ) for our N t = 12 simulations, so we have not ruled out finite volume effects. 
The mass dependence of the position of the peaks in chiral susceptibilities has not been 
adequately explored to check that we are really seeing the chiral limit. We have only used 
one action, and have not explored whether improved actions such as those used by the Lattice 
Higgs Collaboration and Degrand, Shamir and Svetitsky might show different behaviour. 

While the most direct way of answering these questions would be to continue our work to 
lattices with larger N t values and smaller quark masses, such simulations would be expensive, 
and it is not clear if repeating our studies with N t = 16 or 18 would provide the desired 
clarification. There are, however, other less expensive studies which could potentially help 
clarify the situation. The first, which we are pursuing, is to extend our simulations of QCD 
with 3 colour-sextet quarks to N t = 12. This theory is almost certainly conformal, and 
even if it is not, the perturbative evolution of its coupling is extremely slow (this is because 
asymptotic freedom is lost at Nf = for QCD with Nf sextet-quark flavours. Hence we 
should expect to see essentially no change in the value of f3 x between N t — 8 and N t — 12. If 
this is observed it would indicate that the Nf = 2 and Nf = 3 theories are behaving rather 
differently as would be expected if the Nf = 2 theory walks. 

One method of testing how well some of the different choices of improved couplings 
work with sextet quarks and its 2 length scales would be to perform an extensive study of 
the position of the chiral transition with sextet quarks in the quenched theory. That this 
transition is separated from the deconfinement transition was shown in very early studies 


The advantage of this approach is that the production of very large quenched lattices can 
be performed very cheaply, and multimass inversions, already used in the RHMC algorithm 
will allow us to study the chiral condensate over a large range of masses. Here it will be 
possible compare the evolution of [3d and [3 X with various improved couplings. 

Another direction we are pursuing (we thank Julius Kuti for this suggestion) is to simulate 
QCD with 2 colour-sextet quarks at a fixed [3 value above /3 x (N t = 12) - we choose (3 = 6.9, 
- and simulate on lattices with a fixed spatial volume, varying N t to look for the transition. 
Perturbation theory predicts that (3 X = 6.9 for some N t in the range 18 < Nt < 20 . We 
are simulation on 24 3 x N t lattices with 8 < A') < 24 with masses as low as m — 0.00125. 
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Indications are that we will need to increase the spatial box size to accommodate the larger 
N t s, since we are seeing finite volume effects at N t = 22 and 24. 

Temporarily ignoring questions as to whether this theory is QCD-like, we have started 
zero-temperature simulations at /3 = 5.81 (/ 3 d for N t = 12), on 24 1 2 3 4 5 6 7 x 48 lattices. One of the 
reasons for this parameter choice, is to answer another question posed by Julius Kuti, who 
asked what value we estimate for T d / f n . So far we have produced 250 lattices separated by 
100 trajectories for m = 0.01 (m n ~ 0.25) and 250 lattices for m = 0.005 (m T ~ 0.175). 
Larger-lattice zero-temperature simulations at weaker couplings are being contemplated. 

Note that all the simulations reported in this paper are in the state where the argument 
of the Wilson Line (Polyakov Loop) is close to zero. Those states where the argument of the 
Wilson Line is close to ±27 t/ 3 are being ignored, except for /3s approaching (3d and below, 
where transitions between these 3 states are frequent. However, should these be the true 
vacua, charge conjugation would be spontaneously broken ( 9 ] , and so presumably would CP. 
If so, this could possibly provide a mechanism for baryogenesis. 
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